Absolute properties of the low-mass eclipsing binary CM Draconis 
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ABSTRACT 

Spectroscopic and eclipsing binary systems offer the best means for determining accurate 
physical properties of stars, including their masses and radii. The data available for low-mass 
stars have yielded firm evidence that stellar structure models predict smaller radii and higher 
effective temperatures than observed, but the number of systems with detailed analyses is still 
small. In this paper we present a complete reanalysis of one of such eclipsing systems, CM Dra, 
composed of two dM4.5 stars. New and existing light curves as well as a radial velocity curve are 
modeled to measure the physical properties of both components. The masses and radii determined 
for the components of CM Dra are Mi = 0.2310 ± 0.0009 Mq, Ma = 0.2141 ± 0.0010 M©, 
Ri = 0.2534 ± 0.0019 R©, and i?2 = 0.2396 ± 0.0015 R©. With relative uncertainties well below 
the 1% level, these values constitute the most accurate properties to date for fully convective 
stars. This makes CM Dra a valuable benchmark for testing theoretical models. In comparing 
our measurements with theory, we confirm the discrepancies reported previously for other low- 
mass eclipsing binaries. These discrepancies seem likely to be due to the effects of magnetic 
activity. We find that the orbit of this system is slightly eccentric, and we have made use of 
eclipse timings spanning three decades to infer the apsidal motion and other related properties. 

Subject headings: binaries: eclipsing — binaries: spectroscopic — stars: late-type — stars: fundamental 
parameters — stars: individual: CM Dra 
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1. Introduction 

Late-type stars are the most common objects 
in the Galaxy, yet their fundamental properties 
are still not well understood, in part because their 
accurate measurement is challenging. Double- 
lined eclipsing binary systems (hereafter EEs) 



lege, St. Davids, PA 19087, USA 

^°Planetary Science Institute, 1700 E. Fort Lowell Rd., 
Suite 106, Tucson, AZ 85719, USA 

^'^NASA Postdoctoral Program Fellow, Goddard Space 
Flight Center, 8800 Greenbelt Rd Code 690.3, Greenbelt, 
MD 20771, USA 

l^Lowell Observatory, 1400 West Mars Hill Road, 
Flagstaff, AZ 86001, USA 

Department of Astronomy, University of Wisconsin- 
Madison, 475 North Charter Street, Madison, WI 53706, 
USA 



have proven to be the best source of accurate 
properties for fow-mass stars, and a number of 
those sy stems have aheady been studied in de- 
tail (see Ribad 20061 for a review). These anal- 
yses have revealed that low-mass stars in EBs 
have radii that are ~10% larger and effective tem- 
peratures that are ~5% cooler than the predic- 
tions of stellar structure models. On the other 
hand, their luminosities seem to agree well with 
model calculations. These discrepancies have 
been attributed to the effects of magnetic acti-v 
ity o n the component stars (e.g., [Torres & Ribas 



20021: iLopez-Morales &: Ribad l2005t iTorres et al 
20061: iLopez-Moraled l2007t [Morales et al.l l20o'i 

Ribas et al.l 2008l )~ Additional systems with ac- 
curately known stellar properties that cover the 
entire range of sub-solar masses are needed to bet- 
ter constrain the differences between models and 
observations. 

CM Draconis (hereafter CM Dra, GJ 630. lA, 
aj2ooo.o = 161^34-20.^35, <5j2ooo.o = +57°09'44."7) 
is a V = 12.9 mag EB system at a distance of 
14.5 pc from the Sun, which forms a common 
proper motion pair with a,V — 15 mag white dwarf 
(GJ 630.1B, aj2ooo.o = 16h34-21.''57, Jj2ooo.o = 
-h57°10'09."0) at a separation of ~26 arcsec. This 
common proper motion pair moves at a relatively 
large angular speed of roughly 2 arcsec per year, 
which may be indicative of Population II mem- 
bership. Because of this, it has been considered a 
useful system for estimating the primordial helium 
abundance of the Unive rse through th e detailed 
study of its components ( Paczvnski &: Sienkiewic^ 
1984 



CM Dra was first inve stigated spec troscopically 
and photo metrically by Lacy (| l977 ). and more 
recently by iMetcalfe et al. (1996). Both studies 
indicate the system is composed of two similar 
dM4.5 stars with masses of about 0.23 and 0.21 
MfD, orbiting each othe r with a period of 1.27 days. 
Viti et al. I ()1997L I2OO2I) estimated a metallicity of 
-1.0 < [M/H] < -0.6 for the system, and in- 
ferred an effective temperature of 3000 < Toff < 
3200 K. In this paper we describe new observa- 
tions of this binary that add significantly to the 
body of existing measurements. The unique posi- 
tion of CM Dra as the best known binary system 
composed of fully convective stars makes it ex- 
ceptionally important for testing models of such 
objects, and fully justifies a reanalysis in the light 



of our new observations. 

An important feature of this system is that, un- 
like other well-known low-mass EBs, its orbit has 
a small but measurable eccentricity. The precise 
measurement of eclipse timings over a sufficiently 
long period can thus potentially lead to the detec- 
tion of apsidal motion. The rate of this motion can 
be used to infer the internal structure constant k2 
(|Kopallll978Hciaret fc Cimenedl 19931 ). with which 
further tests of the models are possible. In addi- 
tion, the investigation of the times of eclipse can 
also reveal the presence of third bodies in the sys- 
tem through the light-time effect. Attempts to 
detect planets around CM Dra in this way have 
been carried out in the context of the Transit 



of Extrasolar Planets Project fTEP. iDeeg et al . 
1998L I2OOOI : iDovle et al.l I2OOOI : lOeeg et al.l l2008h . 



although no compelling evidence of such objects 
has been found as yet. 

In this paper we present a thorough reanalysis 
of CM Dra to determine the fundamental proper- 
ties of its components, including the masses, radii 
and effective temperatures. Additionally, we have 
measured the rate of advance of the line of ap- 
sides, which turns out to be dominated by the 
General Relativity contribution. In the following 
we describe first all available photometric and ra- 
dial velocity measurements. A combined analysis 
of all the informatio n using th e Wilson-Devinney 
code (hereafter WD, iWilson fc Devinnevlll97l[) is 
discussed in §[31 and in §|3]the times of minimum 
are used to estimate the apsidal motion. The ab- 
solute properties of each component and the age 
and metallicity of the system are derived in §[5l 
and compared with stellar model predictions in 
§[51 Finally, we summarize our conclusions in §[71 

2. Time-series data for CM Dra 
2.1. Light curves 

The photometric data available for CM Dra 
come from a variety of sources. In addition to 
making use of the or i ginal l ight curve in the / band 
published by iLacvl (I1977I) . six new light curves 
measured in the / and R bands have been ob- 
tained with the 0.8m Four College Automatic Pho- 
toelectric Telescope (hereafter FCAPT) located at 
Fairborn Observatory in southern Arizona in the 
Patagonia Mountains. Differential photoelectric 
photometry was conducted from 1995 - 2005 on 
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Table 1: Differential R- and /-band photometry for 
CM Dra from FCAPT. Dates are given in helio- 
centric Julian days on the TT time scale (HJED). 



Table 2: Differential r-band photometry for 
CM Dra from Sleuth. The full version is available 
electronically. 
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available electronically. 
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335 nights. The photometry was typically con- 
ducted using the Cousins R and I filters. The pri- 
mary comparison and check stars were HD 238580 
and HD 238573, respectively. Integration times of 
20-sec were used and the typical precision of the 
delta-R and -I band measures was 0.014 mag and 
0.011 mag, respectively. The relatively large un- 
certainties arise mainly from the faintness of the 
CM Dra and uncertainties in centering the vari- 
able star using blind-offsets (rather than direct ac- 
quires) . 

An additional light curve in the r-band was 
gathered with the Sleuth telescope located at 
the Palomar Observatory in southern Califor- 
nia. Sleuth was one of three telescopes that 
together made up the Trans-atlantic Exoplanet 
Survey (TrES), and its primary use was to dis- 
cover transiting p lanets orbiting stars brighter 



than V = 13 



O'Donovan et al. 2006a. 2007 



Mandushev et al. I I2OO7I) . Sleuth consists of a lens 
with a physical aperture of 10 cm that images 
a field of view of size 5.7 degrees-square onto a 
thinned, back- illuminated CCD with 2048x2048 
pixels, corresponding to a plate scale of 10 arc- 
seconds per pixel. From UT 2004 March 29 to 
UT 2004 June 6, Sleuth observed (as part of its 
survey for transiting planets) a field centered on 
the guide star HD 151613, and this field fortu- 
itously contained our target CM Dra. Whenever 
weather permitted operation, the telescope gath- 
ered exposures in r-band with an exposure time 



of 90 s and a CCD readout time of 27 s, for a 
cadence of 117 s. We used a photometric aperture 
of radius 30 arcseconds (3 pixels) to produce the 
differential photometric time series listed in Table 
m and shown in Fig. H The calibration of TrES 
images, the extraction of the differential photo- 
metric time series (based on image subtraction 
methods), and the decorrelation of the resulting 



light c urves is described e lsewh ere CDunha m et al 
2004 iMandushev et al] l2005t lo^'Donovan et al 



2006bD . The FCAPT and Sleuth data are col- 
lected in Table [Hand Table H 

Collectively these light curves cover the observ- 
ing seasons 1975, 1996-2001, and 2004, and add 
up to more than 20000 individual measurements. 
The short orbital period of CM Dra, along with 
its possibly old age (see below), make it very likely 
that tidal interactions have forced its components 
to rotate synchronously with the orbital motion. 
Thus, although the stars are fully convective, it is 
not surprising that they show a high level of chro- 
mospheric activity as ind icated, for instance, by 
the presence of fl a res (e.g., Eggen fc Sandagell967t 
Lacv et al.l[l976l : iNelson fc Catonll2007^ . Surface 



features (spots) are also conspicuously present and 
are responsible for modulations in the light curves 
that change from season to season. This compli- 
cates the analysis significantly. Prior to combining 
the different data sets, it is therefore necessary to 
correct the light curves for these distortions. Ad- 
ditionally, the large proper motion of CM Dra on 
the sky is such that the system is approaching an 
R = 16.5 mag star to the NW, as shown in Fig. [1] 
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CM Dra 



Fig. 1.— POSS-II DSS2 image in the R band 
showing the position of CM Dra at epoch 1991.5 
and its proper motion on the sky. The common 
proper motion white dwarf companion is labeled 
"WD". Tick marks on the path of CM Dra are 
given in steps of 10 years from 1970 to 2010. North 
is up and East is left. 



Because different photometric apertures have been 
used to obtain the measurements, the proximity 
of this star implies that the light curves from dif- 
ferent instruments may be affected to different de- 
grees by third light. This contamination must also 
be removed before the data can be combined. 

The correction for these spot effects and third 
light contribution was performed by carrying out 
preliminary fits to the light curves in each individ- 
ual data set using the WD code. This program as- 
sumes a relatively simple spot model in which the 
features are circular and uniform. Nevertheless it 
is adequate as a first-order description. FCAPT 
R- and J-band data from the same season were 
used simultaneously. In the absence of spots, the 
model parameters in these fits that depend only 
on light curves are the eccentricity (e), the initial 
argument of the periastron (w), the inclination an- 
gle (i), the temperature ratio (Tcs,2/ToS,i, where 
subindex 2 indicates the less massive component), 
the surface pseudo-potentials (fli), which are re- 
lated to the relative radii (r^), and the passband- 
specific luminosity ratio {j^)- Properties that rely 
on the radial velocities, i.e., the semimajor axis 



(a), the mass ratio (M2/M1) and the systemic 
radial ve locity (7), were held fixed at the values 
given bv iMetcalfe e rZI (Il996h . Limb darkening 
coefficients for these WD runs were computed for 
the standard Cousins R and / bands implemented 
in the code to account for possible corrections of 
these coefficients according to stellar properties at 
each iteration. The stellar atmosphere files in our 
WD implementation do not consider the Sloan r 
band. Therefore, we carried out a number of tests 
to check the adequacy of assuming Cousins R band 
for the Sleuth light curve. We did so by consider- 
ing solutions incorporating the differential in the 
limb darkening coefficients between the Cousins R 
band and the Sloa n r' band (com patible to Sloan 
r) calculated from Claretl (2004). No significant 
effects were found, owing to the fact that the light 
curve shape is quite insensitive to small variations 
in the darkening coefficients in the red bands. 

The spot properties that can be solved for us- 
ing the WD code are their angular radius (rg), 
the longitude {4>s), latitude (Og), and tempera- 
ture contrast relative to the photosphere (Tg/Teg, 
where Tg represents the spot temperature). Be- 
cause 6s, Tg, and Ts/Tcs are strongly correlated 
and can usually not be determined all at once, 
the procedure to fit the spots was carried out by 
iterations. We first computed solutions with vari- 
able (t)s and for several fixed values of , trying 
both dark and bright spots with moderate temper- 
ature ratios. Several scenarios with spots on one 
or both components were tested. Upon reaching 
convergence, we fixed the values of 0s and rg and 
solved for 6s and Ts/T^s- In cases where conver- 
gence was not reached, we selected the fits with 
fixed values of 6s and Ts/T^s yielding the smallest 
residuals. In all of these fits third light (£3) was 
considered as a free parameter as well. The solu- 
tions for the spot parameters and third light that 
give the smallest residuals are shown in Table [3l 
As seen, the spot configurations change somewhat 
from season to season, providing some evidence 
of either redistribution of the features or appear- 
ances and disappearances. Third light is also seen 
to vary from data set to data set for the reasons 
indicated above. In particular, the much larger 
value for the Sleuth data reflects the large pixel 
scale of that instrument, which makes contamina- 
tion by neighboring stars more likely. 

The corrections for spots and third light in each 
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Table 3: Spot and third light parameters from fits to the light curves in each season. Third light is given as 
the percentage of the total light coming from the system at phase 0.25. Parameters labeled as fixed were 
obtained from the trial fits giving the best residuals. 
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Fig. 2. — Differentia l effec t of st ar spots on the 
J-band hght curve of iLacvl (jl977[ ). The sohd Hue 
represents the model described in the text. 



mation resides almost completely in the eclipse 
phases (it depends mostly on their detailed shape) 
so that averaging outside of eclipse has essentially 
no impact on the results. We therefore averaged 
the observations outside of eclipse from the same 
instruments into bins of 0.04 in phase. This pro- 
cedure was applied to the FCAPT observations 
and the Sleuth observations. The total number of 
points used in the solutions is 5356. Unit weight 
was assigned to observations that have no reported 
errors, as is the case for the FCAPT data and also 



Lacvl (|1977l ). whereas individual weights were used 
for the Sleuth observations, for which internal er- 
rors are available. For the out-of-eclipse averages 
from FCAPT and Sleuth we adopted as weights 
the number of combined points and the reciprocal 
of the standard deviation squared, respectively. 



season were computed from the difference between 
the theoretical curves from these full fits and syn- 
thetic curves calculated with the same geometric 
and radiative parameters but with no spots and no 
third light. We then subtracted these effects from 
the original data. As an example, Fig.[2]shows the 
differ e ntial effect of the spots for the light curve of 
Lacvl (Il977h . 

With these transitory effects removed, the pho- 
tometric data can be combined more easily for 
analysis with WD. For practical reasons, we found 
it convenient to bin the large number of original 
data points in order to reduce the computing time 
for the light curve solutions. The relevant infor- 



2.2. Radial velocity data 

For the present study we have made use of 
the sa me spe ctroscopic material discussed by 



Metcalfe et al. 



(fl99(j ). obtained over a period of 
nearly 5 years with an echelle spectrograph on the 
1.5m Tillinghast reflector at the F. L. Whipple 
Observatory (Mount Hopkins, Arizona). These 
observations were taken at a resolving power 
A/AA ~ 35,000, and cover approximately 45 A in 
a single order centered near the Mg I b triplet at 
^^5187 A. Fo r further details we re fer the reader to 
the work of iMetcalfe et al. ( 19961 ). Here we have 



reanalyzed these spectra with improved techniques 
compared to the original study. Radial velocities 
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were obtained with TODCOR ( Zucker fc Mazeh 



19941 ). a two-dimensional cross-correlation algo- 
rithm. The template for both components was 
chosen to be an observation of Barnard's star 
(GJ 699, M4Ve) taken with a similar instrumental 
setup, which provides a close match to the spectral 
type of CM Dra. Unlike the original study, here we 
have made a special effort to match the rotational 
broadening of each component by convolving the 
spectrum of Barnard's star (assumed to have negli- 
gible rotation) with a standard rotational profile. 
The values of the projected rotational velocity 
of the components (vsini) that provide the best 
match to the stars are 9.5 ± 1.0 kms~^ for the 
primary and 10.0 ±1.0 kms^^ for the secondary. 
The average light ratio derived from these spectra 
is L2/L1 = 0.91 ± 0.05 at the mean wavelength of 
our observations. 

As a test, we experimented with other tem- 
plates obtained with the same instrumentation to 
investigate the possibility of systematic errors in 
the velocities due to "t emplate mismatch" (see. 



e.g.. Griffin et al.l l2000[ ). which might bias the 



mass determinations. The use of a template made 
from an observation of the star GJ 725 A (M3.5V) 
produced rather similar velocities, and an orbital 
solution with nearly identical elements and for- 
mal uncertainties only slightly higher than our 
previous fit. The minimum masses from this so- 
lution were smaller than our previous results by 
only 0.23% and 0.14% for the primary and sec- 
ondary, respectively, which are below the formal 
errors in those quantities. A template from an 
observation of GJ 51 (M5.0V) gave an orbital so- 
lution that was significantly worse, and minimum 
masses 0.67% and 0.72% higher than those from 
our reference fit. As a measure of the closeness of 
the match to the real components of CM Dra, we 
computed for each template the cross-correlation 
value from TODCOR averaged over all exposures. 
Both of the alternate templates, which bracket the 
spectral type of CM Dra, gave average correlation 
values that were lower than we obtained with the 
GJ 699 template (particularly for GJ 51), indicat- 
ing the match is not as good. The results using 
Barnard's star are thus preferable, and the above 
tests indicate template mismatch is unlikely to be 
significant. 

The spotted nature of the CM Dra implies 
the possibility of systematic effects on the mea- 



sured radial velocities that could bias the inferred 
masses and radii of the stars. In principle the WD 
code can approximately take into account these ef- 
fects in solutions that use spectroscopic and pho- 
tometric observations simultaneously, as long as 
those observations are contemporaneous. Unfor- 
tunately, this is not the case here, and as seen in 
Table [3] the properties of the spots change signif- 
icantly with time. In order to at least provide an 
estimate of the effect, we have performed exper- 
iments in which we perturbed the individual ve- 
locities by adding the radial-velocity corrections 
that WD computes for each of the spot configura- 
tions in Table [3l We then carried out Keplerian 
fits in each case, and we compared them. The dif- 
ferences in the key parameters (i.e., the minimum 
masses Msin'^i, projected semimajor axis asini, 
e, cj, and M2/M1) are always within the formal er- 
rors. This is not surprising, given that the individ- 
ual velocity corrections are typically smaller than 
0.2 kms^^. Nevertheless, to be conservative, we 
have taken half of the maximum difference in each 
parameter as a measure of the possible systematic 
effect due to spots, and added this contribution in 
quadrature to the uncertainties determined from 
the analysis described in §[31 

The measured radial velocities in the heliocen- 
tric frame are listed in Table SI without any cor- 
rections. They s upersede the m easurements re- 
ported by iMetca lfe ct al.] (|l996l ). The median 
uncertainties are approximately 1.2 kms~^ and 
1.4 km s~^ for the primary a nd secondary, respec- 
tively. IMetcalfe et al.l (jl996l ) did not report indi- 
vidual errors for their radial velocities, but we may 
take the rms residuals from their orbit as repre- 
sentative values. Compared to those (1.77 kms~^ 
and 2.33 km s~^ for the primary and secondary, re- 
spectively), our velocities give significantly smaller 
residuals (1.30 kms'^ and 1.40 kms'^ Table E]). 
We attribute this to our use of templates that bet- 
ter match the rotational broadeni ng of each com - 
ponent (see above), whereas Metc alfe et al. ( 19961 ) 
used an unbroadened template. 

3. Analysis of light and radial velocity 
curves 

Prior to combining these curves, the times of 
observation were transformed to the uniform Ter- 
restrial Time (TT) scale in order to avoid dis- 
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Table 4: Radial velocity measurements for CM Dra in the heliocentric frame. The full version of this table 
is available electronically. 
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continuities resulting from the more than 30 leap 
seconds that have been introduced in the interval 
spanned by the various data sets. For the anal- 
ysis in this section we used the 2005 version of 
the 2003 WD code although updates in this ver- 
sion do not affect the fitting mode used for de- 
tached binaries. The program models proximity 
effects in detail, although they are negligible for a 
well-detached system such as CM Dra. The reflec- 
tion albedos for both components were held fixed 
at the value 0.5, appropriate for convective en- 
velopes, and a gravity br ighteni n g coe fficient of 
0.2 was adopted following [Claretl (|2000l ). For the 
limb darkening we adopted the square root law, 
with coefficients computed dynamically at each 
iteration from the phoen ix atmosphere models 
(lAllard fc Hauschildtl[l995h . in order to follow the 
evolution of the components' properties. 

The light and velocity curves were adjusted si- 
multaneously with WD solving for the epoch of 
primary eclipse (To), the eccentricity (e), the ar- 
gument of the periastron (w), the inclination (i), 
the semimajor axis (a), the systemic radial veloc- 
ity (7), the mass ratio {M2/M1), the secondary 
effective temperature (Toff, 2), the luminosity ra- 
tio at each bandpass {L2/L1), and the surface po- 
tentials (rij). To first order the light curves are 
only sensitive to the temperature ratio of the com- 
ponents. Because the limb darkening coefficients 
need to be interpolated from theoretical tables, we 



assumed Tei 



3 100 K according to the results 



of IViti et all (|l997l . 120021 ) and fitted for the value 

of Teff,2- 

Given that the data span over 30 years, we 



initially attempted also to estimate the period 
(P) as well as the apsidal motion rate (w) with 
WD directly from the light curves, simultaneously 
with the other adjustable quantities. We found 
that this did not yield satisfactory results, and 
the value for to was statistically insignificant com- 
pared to its large error. We then chose to set 
uj to zero and fit each of the light curves sepa- 
rately in order to minimize the effects of possible 
changes in lo from epoch to epoch. The period 
was held fixed at the value found in the analysis 
of eclipse timings described later in §31 which is 
P = 1.268389985 ± 0.000000005 days (this value 
does account for the small effect of u!, as described 
below) . We solved for the parameters of each light 
curve, then computed the weighted averages, and 
subsequently solved for the parameters of the ra- 
dial velocity curves. This was iterated until con- 
vergence, as judged by the changes from one it- 
eration to the next compared to the internal er- 
rors reported by the WD code, i.e., convergence is 
reached when corrections are smaller than errors. 

Table [5] presents the model fits to the differ- 
ent data sets, with the results from all FCAPT 
seasons combined into a single solution. The fi- 
nal column lists our adopted solution in which we 
have taken the weighted average of each parame- 
ter, with weights assigned according to the RMS 
residuals of the fits. The formal errors reported for 
the averages are our more conservative estimates, 
computed as the quadratic sum of the standard 
deviation from the different fits and the internal 
maximum (statistical) error given by the WD code 
for each parameter. The parameters from the fit 
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Table 5: Light-curve solutions for CM Dra from the different data sets. 
1.268389985 days. 



The period adopted is P = 



Parameter 



Lacy 



FCAPT 



Sleuth 



Average 



Physical properties 
To (HJD-2400000) 

■^0 (°) 

in 

Hi 



42958.620510(24) 
0.00521(56) 
108.1(2.2) 
89.784(64) 
15.736(50) 
15.631(59) 
0.06757(12) 
0.06350(17) 



51134.661970(13) 
0.00686(50) 
101.9(0.9) 
89.770(28) 
15.877(39) 
15.506(40) 
0.06700(12) 
0.06403(12) 



53127.302690(21) 
0.00424(56) 
113.9(3.8) 
89.712(62) 
15.862(61) 
15.582(75) 
0.06690(17) 
0.06377(17) 



48042.32743(24)=' 
0.0054(13) 
107.6(6.3) 
89.769(73) 
15.79(11) 
15.59(10) 
0.0673(5) 
0.0637(4) 



Radiative properties 

Teff.J (K) 

(i2/il)H, (ph. 0.25) 
(L2/Li)i (ph. 0.25) 
Albedo 

Gravity darkening 



0.9984(7) 
0.8764(43) 



3100 (fixed) 
0.9926(4) 
0.8721(32) 
0.8782(33) 
0.5 (fixed) 
0.2 (fixed) 



Limb darkening coefficients (square root law) 

XI k yi R 0.268 & 0.690 

X2 Sz yi R 0.293 & 0.669 

XI kyi I -0.043 & 1.011 

X2 k y2 I -0.018 & 0.991 



0.9923(5) 
0.8632(63) 



0.9960(40) 
0.8654(89) 
0.8768(44) 



Other quantities pertaining to the fits 
an (mag) 

aj (mag) 0.0071 
A^obs 830 



0.0236 0.0137 
0.0130 

1656 (R) , 1691 (I) 1179 



"Reference epoch of each light curve corrected to a central epoch. 
''Volume radii. 



to the radial velocity curves are listed in Table [6] 
The results for the eccentricity and co are consis- 
tent with those derived from the light curves. The 
fitted light and velocity curves are shown in Fig. [3] 
and Fig. g) 

The parameters from our light c urve f i ts are 
gener ally similar to t hose reported bv lLacvl (1973) 
(and iMetcalfe etall [1996. who adopted Lacy's 
photometric results), with the exception of the rel- 
ative radius for the secondary, r2. Our value is 



2.1% larger than that determined bv lLacvl (119771 ). 
This discrepancy is significant, corresponding to 2 
times the combined uncertainities. One possible 
explanation is numerical differe nces in the model- 
i ng te chniques: Lacy used the iRussell fc Merrill 
( 19521 ) method, whereas we used WD. Another 
is the treatment of the spots: Lacy assumed the 
spot modulation to be sinusoidal, whereas we per- 
formed a more sophisticated modeling with WD. 



Significant differences in the shape of the mod- 
ulation occur near the eclipse phases, as shown 
in Fig. 121 which can influence the detailed shape 
of the eclipses on which the relative radii de- 
pend. Additionally, Lacy considered the orbit of 
CM Dra to be circular, whereas it is now known 
to be slightly eccentric. Because of the impact 
of r2 on the absolute dimensions of the binary, 
we have investigated this difference by perform- 
ing a number of light-curve solutions based on 
Lacy's data. For this we used the modehng code 
EBOP (Eclipsing Binary O rbit Program; lEtzell 
Il98ll : IPopper fc Eteelll98lh . Under the same as- 
sumptions adopted by Lacy (sinusoidal correction 
for spots, circular orbit) we obtain results very 
close to his for all parameters, indicating the nu- 
merical technique for the modeling is relatively 
unimportant. For a circular orbit but a rectifi- 
cation for spots computed with WD, as we have 
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Phase Phase Phase 

Fig. 3. — Left: Rectified light curves of CM Dra after subtracting the effects of third light and spots, 
separated by instrument. Observations outside of eclipse are binned as described in the text. Residuals are 
shown at the top. Right: Enlargement around the eclipse phases. All light curves are plotted as differential 
magnitude vs. phase, and residuals are plotted in the same order as the light curves. Note the different scales 
for the residuals of each instrument. 




-0.2 -0.1 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Phase 



Fig. 4. — Phase-folded radial velocity curves of 
CM Dra, with the primary shown with filled sym- 
bols and the secondary with open symbols. Resid- 
uals are shown at the top, and the elements of the 
spectroscopic fit are given in Table [6l 

done in our own fits, the results differ somewhat 
from Lacy's, particularly in r2 but also slightly in 
the sum of the relative radii. The largest differ- 
ence, however, is seen when abandoning the as- 
sumption of a circular orbit. We conclude that 
this effect, with some contribution from the treat- 
ment of s pots, has int roduc ed subtle biases in the 



results of iLacvl (|l977[) and iMetcalfe et all (|l996l ) 



that are avoided in the present analysis, and gives 
us confidence in the accuracy of the absolute prop- 



erties described below. 

4. Analysis of the times of minimum 
4.1. Apsidal motion 

As mentioned in the previous section, CM Dra 
has a very small but significant orbital eccentric- 
ity, seen not only in the light curves but also in 
the radial-velocity curves. Both General Relativ- 
ity and the classical theory of tides predict that 
a close system such as this should experience a 
certain degree of periastron advance. Despite our 
attempts described earlier, we were unable to de- 
tect a significant apsidal motion rate (cO) in our 
light curve solutions, even though those data span 
nearly 30 years. However, additional informa- 
tion is available in the form of eclipse timings for 
both minima, and we examine these measurements 
carefully below to investigate possible changes in 
the separation between the primary and secondary 
eclipses that would be indicative of apsidal motion. 

Numerous eclipse timings for CM Dra have 
been reported in the literature usin g a va r iety o f 
techniques, beginning with those of iLacvl (|l977l ). 
Photoelectric or CCD measurements have greater 
precision and are the most useful for our purposes. 
Several timings were obtained in the FCAPT and 
Sleuth observation campaigns and additionally, 
new timing measurements have been made here 
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Table 6: Spectroscopic solution for CM Dra. 
Period is held fixed at the value P = 
1.268389985 days. 



Parameter 



Value 



Physical properties 
To (HJD) 

^0 n 

Ki (kms-i) 
K2 (kms-i) 
a 

7 (kms~^) 
M2/M1 



2446058. 56471±0.00026 
0.0051±0.0013 
129±16 
72.23±0.13 
77.95±0.13 
3.7634±0.0046 
-118.24±0.07'' 
0.9267±0.0023 



RMS residuals from the fits 
Primary (kms~^) 
Secondary (kms^^) 



1.30 
1.40 



"De-projected by adopting the inclination angle from the light 
curve solutions (see Table [Sjl. 
''The true uncertainty of 7 may be larger due to external errors. 



with a number of telescope facilities, as follows. 

A total number ot 20 minima were obtained at 
the Ondfejov observatory with the 65-cm reflect- 
ing telescope with the Apogee AP-7 CCD camera 
in primary focus. The measurements were done 
using the Cousins R filter with 30 s exposure time. 
The nearby star GSC 3881.1146 on the same frame 
was selected as a primary comparison. No correc- 
tion for differential extinction was applied because 
of the proximity of the comparison stars to the 
variable and the resulting negligible differences in 
airmass. The new precise times of minima and 
their errors were determined by fitting the light 
curves with polynomials. 

8 CCD minima were obtained during 2007 and 
2008 in the Sloan r' band using the 2.0m Liver- 
pool Telescope in La Palma. High quality pho- 
tometry (3-4 mmag per image) was obtained, with 
typically 100 photometric points per event. 63 
CCD minima were obtained at the Bradstreet Ob- 
servatory of Eastern University. The equipment 
consisted of a 41-cm f/10 Schmidt-Cassegrain re- 
flector coupled to a Santa Barbara Instruments 
Group ST-8 CCD camera binned so as to give a 
scale of 0.93" pixcl^^. All observations were taken 



through a Cousins / filter. The comparison star 
used was GSC 3881.421 which was always con- 
tained within the same 13' x 9' field. The exposure 
times were 25 sec in duration, typically resulting 
in uncertainties of 3 mmag for each data point. 
Finally, a secondary eclipse of CM Dra was mea- 
sured with the 1.2-m telescope at the F. L. Whip- 
ple Observatory in Arizona using a 4Kx4K CCD 
camera (KeplerCam), binned to provide a scale of 
0.67" pixel" ^. Observations were made through 
a Harris / filter relative to a set of 30 compari- 
son stars, and exposure times were 30 sec each. 
Photometric measurements were performed with 
IRAF using an aperture of 6", and have typical 
uncertainties of 2 mmag. In these three latter 
cas es, times of minima w e re co mputed by using 
the Kwee fc van Woerden ( 19561 ) method. 

All of these measurements (including those 
from the literature) have been converted to the 
uniform TT scale, and are presented in Table [71 
which contains a total of 101 primary timings and 
99 secondary timings. Eclipse timing events com- 
ing from different sources were weight-averaged. 
These data span approximately 35 years, although 
there is an unfortunate gap in the coverage of 
nearly 18 years. 

In the presence of apsidal motio n the times of 

minim um can be described following lGimenez fc Bastero 
(|l995l ) as 



T,=To + P[E 



J-1 



(2j-3)Ai— cost^ + 0(e2). 



(1) 



where j indicates a primary or secondary eclipse 
(1 or 2, respectively), E is the cycle number, and 
Ai is a coefficient dependent on the inclination 
and eccentricity. The first two terms represent the 
linear ephemeris, and the third is the contribution 
of the apsidal motion. Given that the eccentricity 
of CM Dra is very small, powers of or higher 
in these equations have been ignored since they 
produce corrections only of the order of 0.2 sec, 
which are much below the measurement errors of 
the timings. 

For CM Dra we find that Ai k, 2, since the 
inclination is close to 90° and the eccentricity is 
small. Eq. [1] predicts that the deviation of the 
times of minimum from a linear ephemeris has a si- 
nusoidal shape with a semiamplitude of ^188 sec. 
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Table 7: Photoelectric and CCD eclipse timings for CM Dra. The full version of this table is provided in 
electronic form. 



HJED 


(O-C) (s) 


Error (s) 


Prim. /Sec. 


Ref. 


2441855.75476 


-25.4 


30.2 


II 


2 


2442555.90592 


-35.6 


30.2 


II 


2 


2442557.80955 


54.7 


30.2 


I 


2 


2442607.91053 


18.9 


30.2 


II 


2 


2442888.85928 


49.9 


30.2 


I 


2 


2442893.93299 


62.9 


30.2 


I 


2 


2442912.95925 


98.4 


30.2 


I 


2 


2442966.86433 


-30.8 


30.2 


II 


2 


2442994.76890 


-31.6 


30.2 


II 


2 


2449494.63438 


55.1 


2.8 


I 


4 



and a 180° phase difference between the primary 
and secondary. Assuming the rate of periastron 
advance is constant, we may write uj = luo + uj ■ E, 
where E represents the orbital cycle and to is the 
total apsidal motion of the system. The latter 
can be determined from fits of Eq. [T] to each type 
of timing measurement. In the approximation of 
small values of ujE, Eq. [l]can be written as 



T,~To + P[E 
eP 



J-1 



(2) 



+ (2j — 3) Ai — (cos luq — sin luq ■ ujE) . 



A linear fit to the timings can thus be performed 
as 

(O - C)^ = To,, + VjE, (3) 

where Tqj can be taken to represent an effective 
epoch of reference for both minima including the 
effect of the eccentricity, and Pj plays the role of a 
period for each type of minimum. The ephemeris 
for the binary can then be written as 

P1+P2 



P 



To = 



%.l + %.2 P 



(4) 
(5) 



Given values for the orbital elements {P, e) and 
Ai , we may compute luq and uj from the linear fit 
parameters Toj and Vj as 



arccos 



2n %,,-To,,-i_\ 



AieP 2 

2tt Vi - V2 
AieP sinujQ 2 



(7) 



From Eq. [7] it can be seen that a difference in the 
periods for each type of minimum is an indication 
of the presence of apsidal motion in the binary. 

Fig. O shows the O — C values for the primary 
and secondary minima of CM Dra as a function 
of the cycle number. There would seem to be a 
linear trend although the scatter of the measure- 
ments is fairly large. This scatter may be due 
in part to random errors, but there could also be 
biases arising from the presence of spots on the 
surface of the components. As a test, we sim- 
ulated light curves for CM Dra with the different 
spot configurations given in Table [3l and we found 
that the presence of spots can indeed skew eclipse 
timing determinations by up to ~ 1 5 sec. Simi- 
lar results were foimd in a study bv iHargis et al. 
(2000). Because of this effect, observational er- 
rors in the timings could be significantly underes- 
timated. We therefore performed linear fits of the 
times of eclipse with the internal errors doubled, 
in order to preserve the relative weights between 
the measures and obtain a reduced value closer 
to unity. This yielded more realistic uncertain- 
ties for the parameters of the fit. The results are 
shown in Table [8l and indicate an apsidal motion 
of cO = (2.3 ± 1.4) • 10-4 deg cycle'^ i.e., a de- 
tection with 1.6(7 significance. Tests in which the 
internal errors were augmented by adding 15 sec 
in quadrature (to account for the potential effects 
of spots) instead of doubling them gave the same 
results, within the errors. 
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Fig. 5. — Observed minus calculated (O—C) resid- 
uals from the eclipse timings of CM Dra (filled 
symbols for the primary, open symbols for the 
secondary) with respect to a linear ephemeris. 
The linear fits to apsidal motion (solid line) and 
the theoretically predicted apsidal motion (dashed 
line) are shown. 

4.2. Third body effects on the eclipse tim- 
ings 

The analysis of times of minimum can also re- 
veal the presence of third bodies in eclipsing sys- 
tems through the time-delay effect caused by the 
orbit of the binary around the barycenter of the 
system. This produces a sinusoidal m odulation of 
the (O — C) values from the timings. iDeeg et "aT 



(|2008[ ) have recently reported the possible pres- 
ence of a third body around CM Dra based on 
a parabolic fit to their sample of (O — C) values. 
We find, however, that using our own timings a 
parabolic fit is essentially indistinguishable from 
a linear fit to the measurements. Thus, any third 
body must have a period longer than roughly twice 
the span of the measurements, or ~60 years, or 
must induce a light-time effect below '-^ 15 s which 
would be undistinguishable from the dispersion of 
the data due to spot effects. 

Another indication of the possible presence of 
a third body is the small eccentricity of the close 
binary orbit of CM Dra. Systems with periods as 
short as that of CM Dra are usually assumed to be 
tidally circularized early on ((Mazeh 2008). possi- 
bly even during the pre-main sequence phase. To 
explain the present non-zero eccentricity one may 
invoke the presence of a perturbing component in 



Table 8: Results of the linear fits to the eclipse 
timings for apsidal motion. 



Properties 



Weighted fit 



^0,1 (s) 


48042. 32778±0. 00002 


n / 1 — 1 \ 

Vi (s cycle ) 


1.2683899936±0. 0000000064 


xl 


1.303 


%..2 (S) 


48042. 96084±0.00002 


V2 (s cycle^^) 


1.2683899765±0.0000000069 


xl 


0.920 


P (days) 


1.268389985±0.000000005 


To (HJED) 


48042.327214±0.000014 


Wo (deg) 


104.9±3.7 


Lo (deg cycle" ^) 


(2.3±1.4) lO-'^ 


U (years) 


5400±3200 



a more distant orbit. Such a configuration can pro- 
duce secular variations of the orbital parameters of 
the inner orbit, such as an eccentricity modulation 
with a typical period Umod given by 



Umod — P1.2 

ai2 



as \ Ml + M2 



(8) 



where Pi 2 and ai 2 are the period and semima- 
jor axis of the inner orbit of CM Dra, and 03 
and M3 are the semimajor axis of the third body 
around the center of mass of the triple system 
and the mass of the third body, respectively. A 
third body is actually known in the CM Dra sys- 
tem (the common proper motion white dwarf com- 
panion). Adopting a mass for the white dwarf 



of 0.63 M0 from [Bergeron et al.l (|2001r ). along 



with an angular separation from CM Dra of about 
26" (corresponding to ~380 AU at the distance 
of CM Dra), the modulation period on CM Dra 
would be roughly 2 Gyr. However, the effect of 
such a long-period eccentricity pumping would be 
averaged out over many apsidal motion cycles, and 
therefore the orbit would remain circular. One 
may assume that eccentricity pumping by some 
other body in the system will only be effective if 
Umod S 5400 years, which is the period of the ap- 
sidal motion found for CM Dra. This provides 
a constraint on the properties of this putative 
body, if it is to explain the measured eccentric- 
ity. Fig. [5] represents the allowed region (mass vs. 
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Fig. 6. — Allowed region (shaded area) in a mass 
vs. orbital period diagram for a third body in 
the CM Dra system according to observational 
constraints: eccentricity modulation (solid line), 
light-time effect (dashed line) and stability crite- 
ria (dot-dashed line). 



period) of the companion by accounting for the 
non-detection of light-time effect above 15 sec and 
the eccentricity pumping. Also, we consider that 
P^/PcM Dra ^ 30 for Stability criteria of hierar- 
chical triple systems. As can be seen, a massive 
planet or light brown dwarf with an orbital period 
of 50-200 days would fulfill all constraints. 

5. Absolute properties of the components 
of CM Dra 

Based on the fits to the light curves and the ra- 
dial velocities, the absolute physical properties of 
the components of CM Dra including the masses 
and radii can be derived independently of dis- 
tance or flux calibrations. We report these val- 
ues in Table [9] With the measured radii, we find 
that the predicted rotational velocities of the pri- 
mary and secondary, assuming synchronous ro- 
tation, are 10.22 ± 0.08 kms^^ and 9.67 ± 0.07 
kms^^, respectively. These values are in good 
agreement with the v sin i measurements from our 
spectra (§[2l2|). 

The effective temperatures of the components 
are not directly accessible from the light curve 
analysis, which yields only their ratio as measured 
by the relative depths of the eclipses. In §|3] we 
adopted a val ue for the primary Tpff from an ex- 
ternal source (jViti et al.lll997l 120021 ) . based on a 



modeling of the spectrum of CM Dra. It is pos- 
sible, however, to determine the individual tem- 
peratures in another way, using information from 
the light curves along with a combined near-IR 
magnitude for the system and its trigonometric 
parallax. Here we have used the 2MASS mag- 
nitude Ks — 7.796 ± 0.021, subsequ ently con- 
verted to the Johnson system following [Cargenter 



(120011). and the p a rallax tt ~ 69.2 ± 2.5 mas from 
van Altena et al.l ( 19951 ). We chose to rely on 



a near-IR magnitude because the corresponding 
bolometric corrections are less dependent on T^s 
and chemical composition. We began by adopt- 
ing the value of Teg = 3100 K as a starting point 
for the primary, from which the secondary Toff fol- 
lows from the measured temperature ratio. Bolo- 
metric corrections for each star were taken from 
iB essell et al. ( 19981 ) as a function of temperature, 
and averaged since they are virtually identical. 
The total luminosity was then computed. The 
ratio of the luminosities can be calculated from 
the temperature ratio and radius ratio, both of 
which are measured directly and accurately from 
the light curves: 



hi 

Li 



, 2 /TH \ 4 

r-2 \ / Jfi 



ri 



[cff,2 



cff,l 



= 0. 



±0.022. 



(9) 



Individual bolometric luminosities are thus eas- 
ily derived, and since the absolute radii are also 
known, the individual temperatures can be ob- 
tained. This process was iterated until the cor- 
rections to the temperatures were below 1 K. 
The result is independent of the starting point 
for the primary temperature. The mean bolo- 
metric correction resulting from the calculation is 
BCk = 2.66 ± 0.05, and the total luminosity is 
0.0104 ± 0.0009 Lq. The individual temperatures 
and luminosities are listed in Table [9l in which 
the uncertainties include all measurement errors 
as well as an assumed uncertainty of 0.05 mag for 
the bolometric corrections, but exclude systemat- 
ics that are difficult to quantify. The Tefi values, 
which have a m ean of 3125 K, agree ve ry well with 
the estimate of IViti et all (|l997l [20o3 ). 



As a check, we used the above temperatures and 
our light ratios in the R and / bands from the light 
curves to predict the light ratio in V, appropriately 
scaling the NextCen models (jHauschildt et al. 
1999*). The resuh is L2/L1 = 0.86 ± 0.15, which 
is consistent with the value determined spectro- 



13 



Table 9: Absolute physical properties of CM Dra. 
Properties Component 1 Component 2 



M (Mq) 
R (Ro) 
log 5 (cgs) 

Teff (K) 

log(L/LQ) 
Age (Gyr) 
[M/H] 



0.2310±0.0009 
0.2534±0.0019 
4.994±0.007 

3130±70 
-2.258±0.038 



0.2141±0.0010 
0.2396±0.0015 
5.009±0.006 

3120±70 
-2.313±0.056 



4.1 ± 0.8 (Main Sequence) 
-1 < [M/H] < -0.6 



M, 



=4.74 is used to compute luminosities llBessell et al.|[l"998l ). 



Bol©- 



scopically (§ 12. 2|) . within the errors. The mean 
temperature of the system may also be estimated 
from available color indices for CM Dra, and the 
rec ent color /tempera t ure ca libration for M dwarfs 
by ICasagrande et all (l2008l) . We used the Fi? J 
magnitudes of Lacv (SlJ^ and the JHKs mag- 
nitudes from 2MASS to construct twelve different 
color indices (after conversion of the RI magni- 
tudes from th e Johnson syste m to the Cousins sys- 
tem following iLeggett 1992 ) , which are of course 
not independent of each other although they do 
serve to gain a better idea of the scatter among 
the various calibrations. We obtain a weighted 
average temperature of 3050 ± 50 K, which is only 
slightly lower than the estimates above, but has 
the virtue of being completely independent of the 
parallax and the light curve parameters. 

The age, along with mass and chemical compo- 
sition, is an indicator of the evolutionary status 
of a star. When known, it becomes a powerful 
constraint that can be used in the model com- 
parisons. For CM Dra we may obtain a rough 
estimate of its age by considering the proper- 
ties of its whit e dwa, rf companion. According to 
Bergeron et al. I (l200ll) . the cooling age of the white 
dwarf is 2.84 ±0.37 Gyr. Given its estimated mass 
( ^0.63 Mf7)) and t he ini tial-final mass relationship 
of Catalan et al. ( 20081 ) . the mass of the main se- 
quence progenitor is estimated to be 2.1 ±0.4 Mq. 
For stars of such mass, stellar evolution models 



predi ct a hfetime of about 1.3 Gyr (jGirardi et al 
2000l ). We therefore infer an approximate age for 



CM Dra of 4.1 Gyr with a 20% uncertainty level 
coming from uncertainties in the mass of the white 
dwarf progenitor and its metallicity. This total 
age indicates that CM Dra is well on the main 



sequence. 

The chemical composition of CM Dra has been 
notoriously difficult to determine, which is unfor- 
tunate for such an important system. It has usu- 
ally been considered to be metal-poor, although 
this is based mostly on circumstantial evidence 
(i.e., its large space motion). Attempts to deter- 
mine the metallicity by various means have often 
produced inconsiste nt r esults. ICizid (|l997t ) and 
Leggett et al. I (ll998l) concluded the composition is 
near solar, while IViti et al.l (|l997l . |2002[ ) found a 
metal-poor composition (—1.0 < [M/H] < —0.6) 
by performing fits to the spectral energy distribu- 
tion and several diagnostic spectral features us- 
ing stellar atmosphere models. However, some 
systematic differences between the estimates from 
optical and near-IR spectra in the latter studies 
are disconcerting and cast some doubts on the re- 
sults. Our own checks using the same spectro- 
scopic material and the most recent version of the 
NextGen models did not yield an improvement 
in the results. The various metallicity indicators 
still show disagreements, and would seem to indi- 
cate shortcomings in the model atmosphere calcu- 
lations. Thus, the metallicity of CM Dra remains 
poorly determined. 

Kinematics of the common proper motion 
group of CM Dra could provide further insight on 
its age and metallicity. The space velocity com- 
ponents of the system are U = —106.8 kms~^, 
V = -119.8 kms-^and W = -35.1 kms^^. 
These values indicate that the system probably 
does not belong to the thin disk population. No 
clear correlation between kinematics and metal- 
licity or age ha s been found for stars o n the solar 

I2OO4I ). and cer- 



neighborhood ([Nordstrom et al 



tainly no claims can be made on an individual 
star basis. We must conclude that the kinematics 
of CM Dra do not seem to stand in contradiction 
with an age of about 4 Gyr, neither help to discern 
between a solar or moderately sub-solar metallic- 
ity. Another interesting trait of CM Dra is the 
fact that it has remained weakly bound to its rel- 
tively distant companion for a long time. The 
binding energy of the system is over four orders 
of magnitude smaller than its kinetic energy with 
respect to the local standard of rest. Whatever 
perturbations the system has suffered during its 
life, they must come from a smooth potential or 
else the pair would have been broken. This should 
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provide interesting constraints to the mechanisms 
of star acceleration in the Galaxy. 

6. Comparison with theoretical models 

Our mass, radius, and temperature determina- 
tions for CM Dra are compared in Fig. [7] with 
the predictio n s of s tellar evolution models from 



Baraffe et al 



( 19981 ). The measurements support 
a trend found previously for other low-mass EBs, 
in the sense that the observed radii for both com- 
ponents are larger than predicted by theory, in 
this case by ~4.7% and ~5.0% for the primary 
and secondary, respectively. The effective tem- 
peratures are cooler than the models indicate, by 
^^6.8% and ^6.3%. We note also that while these 
are significant offsets (compared to the errors), the 
slope of the models appears substantially correct. 
These deviations refer strictly to the comparison 
with solar-metallicity models, and would be even 
larger if a lower metallicity were assumed. For ex- 
ample, the offsets would increase to ^--^10% in the 
radii for models with [M/H] = — 1. 

Magnetic activity on the components of simi- 
lar low-mass EE systems has often been proposed 
as an explanation for these discrepancies between 
models and observations (see §[T|) . The activity in 
these typically short-period binaries is associated 
with the very rapid rotation resulting from tidal 
synchronization with the orbital motion. One 
manifestation of this activity is the presence of 
surface features (spots) that tend to block a frac- 
tion of the outgoing radiation. The star adjusts by 
increasing its size in order to conserve flux, and at 
the same time the effective temperature becomes 
lower than in a spot-free star. Recent work has 
shown that the same hypothesis appears to ex- 
plain the differences observed between act ive and 
inactive single stars ([Morales et al. I l2008f) . The- 
oretical efforts have had some success in repro- 
ducing the observations for sub-solar mass binary 
systems by accounting for stella r activity in the 
models, at least to first order (e.g.jD'Antona et al 



2000 : iMullan fc MacDonaldl 1200X1: IChabrier et al 
20071 ). In the most recent of these studies the au- 
thors examined the effects of activity in reducing 
the convection efficiency as well as in obstruct- 
ing radiation due to the presence of dark surface 
features. The first of these effects is equivalent 
to a reduction in the mixing length parameter 



(oml), whereas the second can be parametrized in 
terms of the fractional spot coverage. The results 
show that for stars in the fully convective regime 
(M < 0.35 Mq) the effect of a reduction in auL 
is minimal, while the presence of spots has a sig- 
nificant effect, and accounting for this can in fact 
reproduce the properties of CM Dra with a spot 
fraction of about 30%. For somewhat more mas- 
sive stars theory predicts that both a reduced oml 
and spot coverage lead to similar effects on the 
global properties. Although further observational 
and theoretical work is needed, these predictions 
appear at least qualitatively consistent with the 
findings of previous EB studies that suggest the 
radius discrepancies with the models are roughly 
5% for stars with M < 0.35 Mq, and about 10% 
for higher mass stars with convective envelopes. 
The implication is that in the former case the de- 
viations are due only to spots, whereas for stars 
with radiative cores both spots and the reduction 
in the convection efficiency are important. A more 
detailed study of the relationship between activ- 
ity and the radius discrepancies is underway by 
a subset of the present authors, including consis- 
tency checks with all observational constraints for 
late- type stars in binaries. This work will be pre- 
sented in a forthcoming paper, placing CM Dra in 
context with the rest of the low-mass EBs. 

The value of the apsidal motion found in §|4] 
provides a different type of test of internal struc- 
ture models since the rate of classical precession 
induced by tidal effects depends on the density 
profile of the stars. Following the prescriptions in 
( Kopall 19781 ) ■ tidal apsidal motion is given by: 



i^tidal = 360° (C2,lfc2,l + 02,2^2,2) : 



(10) 



where C2a are coefficients that depend on the prop- 
erties of each component, and are the in- 
ternal structure constants dependent on the den- 
sity profile. Using the internal profile models of 



Baraffe et al 



(|1998[ ). we derive logfc2,i = -0.95 
and logfc2,2 = —0.96 for the primary and sec- 
ondary components, respectively. These values 
yield utidai = (1.64 ± 0.04) • 10-^ deg cycle"!. 

However, several more phenomena can con- 
tribute to the magnitude of the total apsidal mo- 
tion, aside from the classical effects of tidal inter- 
action. One is relativistic precession, which de- 
pends essentially on the masses and orbital period 
of the binary system. Following the formula given 
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Fig. 7. — Comparison of the physical properties for CM Dra with models of iBarafFe et al.l (|1998l l in the 

M-R diagram (left) and the M-Tcff diagram (right). Models for different metallicities and ages are shown, 
as labeled. 



by ICimenej (1985), and based on the properties 
derived here for CM Dra, the predicted effect is 
Wrci = (2.711 ±0.005) • 10-4 deg cycle-^ Combin- 
ing the tidal and relativistic contributions, we ob- 
tain a theoretical value of the total apsidal motion 
of Lotkeo = (1-91 ± 0.04) • IQ-^ deg cycle'^ This 
value is incompatible with the observed cuobs = 
(2.3 ± 1.4) ■ 10"^ deg cycle"^. The discrepancy is 
significant, at the 12-cr level. 

A third body may also alter the apsidal motion 
rate of the binary. Interaction with the distant 
white dwarf companion needs in principle some 
consideration. The typical period of the apsi- 
dal motion resulting from a third-body perturba- 
tion is given by Eq. [S] As described earlier, for 
the white dwarf CM Dra companion this modu- 
lation period is about 2 Gyr, thus the contribu- 
tion of the white dwarf to the apsidal motion is 
completely negligible compared to the relativistic 
contribution, which has a much shorter period of 
~4600 yr. There have been some claims of de- 
tections of low-mass companions orbiti ng CM Dra 



from analyses of the eclipse timings (jDeeg et al 
[2000, 2008, e.g.,), but the evidence so far does not 
seem compelling. As we discuss in ^A.'2\ a third 
body may help to explain the small but signifin- 
cant eccentricity and remain undetected to via the 
light-time effect. Given the possible range in mass 
and orbital period of such putative substellar ob- 
ject, we estimate that its contribution to the apsi- 
dal motion of CM Dra may be sufficient (depend- 
ing on the relative orbital inclination) to explain 



the observed differences. Note that scenario has 
been advocated to resolve th e discrepancy found 
in the eclipsing binary DI Her (jCuinan &: Malonev 



19851) 



CM Dra has often been regarded in the past 
as a favorable system for inferring the primordial 
helium abun dance, assuming that it is a Pop ula- 
tion II star. Paczvnski &: Sienkiewicj (1984) de- 
scribed a method using polytropic stellar mod - 
els that was followed by iMetcalfe et al. ( 19961 ). 
who obtained bulk helium abundances of about 
0.3 for both stars. These values are significantly 
higher than estim ates using other methods. An- 
other analysis by IChabrier fc Baraffd (|l995l) led 
to a much lower value of 0.25, through a compari- 
son between their models and the physical pr oper- 
ties of the CM Dra components reported by iLacv 
(|l977l ). Those authors indicated, however, that for 
masses as low as those of CM Dra the models de- 
pend only weakly on the helium abundance. The 
present work shows that standard models are as 
yet unable to reproduce the observed values of R 
and Toff for the components of this binary at their 
measured masses. It would seem, therefore, that 
a determination of the helium abundance by com- 
parison with model predictions (or simpler poly- 
tropes) is not particularly meaningful at the mo- 
ment, given that the significant effects of activity 
on the global properties of these stars are not yet 
properly accounted for. 
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7. Conclusions 

Prompted by significant improvements in eclips- 
ing binary analysis methods since the most re - 
cent major photometric study by iLacvl (Il977h . 
and aided by new photometric observations gath- 
ered here to complement existing data, as well 
as improvements in the spectroscopy, we have 
conducted a thorough reanalysis of the classical 
low-mass double-lined eclipsing binary CM Dra. 
The goal has been to provide the best possible de- 
terminations of the physical properties to enable 
stringent tests of stellar theory. Our results for 
the masses and radii of the stars, which we esti- 
mate to be about 4.1 Gyr old, are Mi = 0.2310 ± 
0.0009 Mq, M2 = 0.2141 ± 0.0010 Mq, Ri = 
0.2534 ± 0.0019 R©, and R2 = 0.2396 ±0.0015 R©, 
with formal relative uncertainties of only ~0.5%. 
A special effort has been made in this study to 
investigate possible sources of systematic error in 
these quantities and to assess their importance. 
We have performed a number of tests during the 
light-curve analysis, the spectroscopic analysis, 
and the determination of effective temperatures. 
The resulting uncertainties of these physical prop- 
erties are thus believed to be realistic, and to offer 
the best opportunity for carrying out meaning- 
ful tests of models of stellar evolution and stellar 
structure for fully convective main-sequence stars. 

We find that the radii and temperatures of 
CM Dra show the same sort of discrepancy with 
model predictions as found previously for other 
low-mass EEs, which are at the level of '~5-7% 
in this particular case. Mounting evidence in- 
dicates that such differences can be ascribed to 
magnetic activity effects. Further research is un- 
derway to estimate the corrections needed in the 
models in order to reproduce the observations of 
low-mass EEs, g iven t he prescriptions proposed by 



Chabrier et al 



( 20071) . 

Measurements of the times of minimum for 
CM Dra clearly show the presence of apsidal mo- 
tion in the system. However, its value is still 
poorly determined on account of observational er- 
rors, other errors due to distortions caused by 
spots, and the Hmited time coverage of the data. 
There also seems to exist a discrepancy bewteen 
the observational value and that derived from 
General Relativity and tidal theory. A third body 
in the system, which may be responsible for the 



non-zero eccentricity, could provide an explana- 
tion to the observed difference. Further measure- 
ments over the coming years will greatly help to 
constrain the precession of the line of apsides and 
separate the two effects more clearly, in addition 
to providing a better basis for investigating the 
possible presence of a third body in the system. 
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